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ABSTRACT 

We report high resolution imaging of the ultraluminous infrared galaxy Arp 
220 at 1.1, 1.6, and 2.22 pm with NICMOS on the HST. The diffraction-limited 
images at 0.1-0.2" resolution clearly resolve both nuclei of the merging galaxy 
system and reveal for the first time a number of luminous star clusters in the 
circumnuclear envelope. The morphologies of both nuclei are strongly affected 
by dust obscuration, even at 2.2 jum : the primary nucleus (west) presents a 
crescent shape, concave to the south and the secondary (eastern) nucleus is 
bifurcated by a dust lane with the southern component being very reddened. 
In the western nucleus, the morphology of the 2.2 ^m emission is most likely 
the result of obscuration by an opaque disk embedded within the nuclear star 
cluster. The morphology of the central starburst-cluster in the western nucleus 
is consistent with either a circumnuclear ring of star formation or a spherical 
cluster with the bottom half obscured by the embedded dust disk. Comparison 
of cm-wave radio continuum maps with the near-infrared images suggests that 
the radio nuclei lie in the dust disk on the west and near the highly reddened 
southern component of the eastern complex. The radio nuclei are separated by 
0.98" (corresponding to 364 pc at 77 Mpc) and the half-widths of the infrared 
nuclei are ~0.2-0.5" . At least 8, unresolved infrared sources - probably globular 
clusters - are also seen in the circumnuclear envelope at radii 2-7" . Their 
near-infrared colors do not significantly constrain their ages. 

Subject headings: galaxies: individual (Arp 220) - Galaxies: active 
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1. Introduction 

Arp 220 (IC 4553/4), with an infrared luminosity of 1.5 x 10 12 L Q at A=8-1000 fan, 
is one of the nearest ultraluminous infrared galaxies ( goiter et al. 1987|) . Visual wavelength 
images reveal two faint tidal tails, indicating a recent tidal interaction (cf. |Joseph fc Wright 



1985|) , and high resolution ground-based radio and near-infrared imaging show a double 



nucleus ( [Baan et al. 19871 , IGraham et al. 1990| ). The radio nuclei are separated by 0.98" at 



PA ~ 90° ( [Baan fc Haschick 1995| ). The system is also extraordinarily rich in molecular 
gas (Mff 2 ~ 9 x 10 9 M & , Scoville, Yun & Bryant 1997) and about 2/3 of this gas and 
dust is concentrated in a thin disk with radii < 250 pc. Ground-based observations at 10 
to 30 /im suggest that the far-infrared luminosity originates from a region of similar size 
( |Wynn- Williams fc Bccklin 1993|) . To power the energy output seen in the infrared by 



young stars requires a star formation rate ~ 10 2 M Q yr -1 . Alternatively, if the luminosity 
originates from an AGN, this source must be sufficiently obscured by dust that even the 
mid-infrared emission lines are highly extinguished since spectroscopy with ISO shows no 
evidence of very high ionization lines at wavelengths out to 40/im ( Sturm et al. 1996|) . 



In this letter we report near infrared imaging at A=l-2.2^m with 0.1-0.2" resolution 
obtained using the NICMOS camera on the HST. These high resolution data clearly reveal 
the complex nature of the nuclei with the morphologies of both near-infrared nuclei strongly 
affected by dust obscuration and enable an improved registration of the infrared relative to 
the radio nuclei. Throughout this letter, we adopt a distance of 77 Mpc (1" corresponds to 
373 pc) for Arp 220 using H =75 km s^Mpc" 1 QSoifer et ai.T987j ). 
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2. Observations and Results 

Arp 220 was observed using the NICMOS camera 2 on the HST on 4 April 1997 
as part of the Early Release Observations (ERO) program. Camera 2 has a 256x256 
HgCdTe array with 0.076" pixels, providing a 19.2" field of view. Images were obtained 
in the F110W, F160W and F222M filters, providing resolutions of 0.11, 0.16 and 0.22" 
respectively. For both the galaxy and an HST guide star, observed to determine the 
point-spread-function (PSF), a four step dither pattern was executed with offsets of 5.4 
pixels (0.41"). At each position in the dither pattern on Arp 220, two integrations, each 
lasting 128 sec, with multiple non-destructive reads (MULTIACCUM) were obtained. The 
data were dark subtracted, flat fielded and corrected for cosmic rays using the IRAF 
software NICMOS-pipeline routines. The dithered images were shifted and coadded using 
the drizzle routine in IRAF. The flux scale employed calibration factors of 2.35 x, 2.77x, 
and 6.11xl0~ 6 Jy (ADU/sec)" 1 at 1.1, 1.6, and 2.22 /an ( Rieke et aJ.T997|) . The measured 



fluxes at H (1.6/mi) and K (2.22/zm) for the nuclear sources and for a 5" aperture agree to 
better than 10% with those obtained from ground-based imaging by Carico et ad. (1990). 
At 1.1 //m a direct comparison is not possible since the standard J filter is at 1.25/mi. The 
rms noise levels in the final images are 12, 15, 35 //Jy (arcsec) -2 at 1.1, 1.6, and 2.22 /mi, 
respectively. 

For observations with NICMOS camera 2, the angular resolutions are a factor of 
two lower at 2.22//m than at 1.1/mi. This can seriously hinder attempts to analyze color 
gradients where there are compact emission sources unless the images at each wavelength 
are convolved or deconvolved to a common resolution. We have adopted two approaches : 
convolving the two short wavelength images to the resolution of the raw image at 2.22/xm 
and deconvolving the two long wavelength images to the resolution at 1.1/mi. The former 
is, of course, the safest approach but provides more limited spatial resolution. For the latter 
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approach, deconvolution was done using the Richardson-Lucy algorithm without damping 
( |Busko fc Stobie 1997|) to achieve 0.14" resolution. 



In Figure 1 (color plate) the images for the entire camera 2 field from the 3 bands are 
combined in false color. The image clearly reveals the structure around the twin nuclei 
separated by ~ 1" E-W in addition to extended emission from the system out to >8". At 
least eight unresolved (< 30pc) sources are also seen at radii 2-7" - the brightest being 
approximately 5" southwest of the nucleus. All of them have bluer colors than the nuclei 
but their magnitudes are sufficiently bright that they cannot be individual stars. We 
hereafter refer to them as clusters. The brightest of these clusters (1, 2, and 4) are seen 
in the visible band HST imaging of Shaya et al. (1994) and they are also clearly seen in 
subsequent I-band HST imaging (Borne and Lucas 1997). 

The central region is shown in contour form in Figures 2-4 (Plates 2-4) using : the 
original drizzled images at each wavelength; images convolved to 0.22" resolution; and 
deconvolved to 0.11" resolution, respectively. In each case, the lower right panel shows 
the ratio of the 2.22/xm/l.l/um data. Coordinate offsets in Figures 2-4 are measured from 
the peak of the deconvolved 2.22/xm emission in the west nucleus. Flux measurements for 
the nuclear components and all identified cluster sources are listed in Table 1 along with 
measured sizes (FWHM) for resolved structures. 



3. Analysis and Discussion 

Most of the near infrared emission from the nuclear region of Arp220 must be stellar 
(rather than hot dust or shock-excited emission lines) since the strength of the CO overtone 
absorption at 2.35/im requires that >90% of the observed near-infrared flux in a 5" aperture 
arises from stars ( Armus et al. 1995 , Shier, Rickc fc Rickc 1996 ). In previous ground-based 
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images (eg. |Graham et al. 1990]) , the core of Arp 220 appears as two nuclei separated E-W 
by 0.95" . In the HST-NICMOS images, the eastern nucleus breaks up into northern (NE) 
and southern (SE) components separated by 0.4" - the latter appears highly reddened and 
is barely visible in the l.l//m image. And, the bright nucleus on the west (W), is seen 
as a crescent arc in the new images. Thus, for both nuclei, their morphology, even at 
A=2.2//m (where the extinction is 10 times less than Ay) appears strongly influenced by 
dust extinction. In addition, a new component is seen ~ 0.6" south of the western nucleus 
(S) and another midway between the two nuclei. 



3.1. Registration 

The accuracy of the guide star position for the Arp 220 observations was insufficient 
to independently determine the absolute positions of the near infrared images. We have 
therefore registered the infrared images by comparison with cm-wave radio continuum maps 
for which the coordinates are well-determined. Graham et al. (1990) matched the brightest 
radio sources (A and B, Paan fc Haschick 1995|) with the brightest portions of the eastern 



and western nuclei (i.e. components NE and W). However, this placement now seems 
unlikely since the separation and position angles of the two radio components and the two 
brightest infrared components are significantly different. Instead, we suggest that both of 
the two brightest radio components are located within the areas of very high extinction - 
south of the W component and between the NE and SE sources. This registration places 
the weak radio source C precisely on the southern (S) infrared source which is relatively 
blue and therefore unlikely to be strongly extincted. Using the 2.22/im/l.l/xm ratio image 
(Figure 3) to estimate the relative extinction across the region, we have crudely corrected 
the 2.22/iin image for extinction under the assumption that the dust is in the foreground, 
i.e. not mixed with the stars. In Figure 5, the radio source contours from Baan & 
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Haschick (1995) are shown superposed on this 'extinction-corrected' 2.22/im image using the 
registration discussed above. Radio source B lies within the dark area below the crescent on 
the W source while A is at the north of the reddened SE infrared source. Support for our 
registration of the near infrared emission is provided by the strong morphological similarity 
of the infrared structure shown here and the H 2 CO emission ( Paan fc Haschick 1995|) - the 
latter exhibits similar offsets relative to the radio continuum. 

3.2. Reddening 

The 2.22//m/l.l//m ratio images in Figures 2-4 may be used to probe variations in the 
extinction assuming that the intrinsic color of the underlying stars is not strongly varying. 
The false color image for the deconvolved data (Plate 1) clearly shows an overall increase 
in the reddening towards the central regions of Arp 220 - in the annulus at 2-4" radii, 
mi. i_2. 2=2. 30 mag, while inside 2" radius , mi.i_ 2 . 2 =3.19 mag. In the central 2" , Figure 3 
reveals even larger variations : high reddening is seen over the SE source where the color 
(1.1-2.2/im) is 3.86 mag and in a nearly complete ring surrounding the bright W source 
where the color is 3.31 mag. For a population of stars with ages < 10 7 yrs, typical colors 
are J-H = -0.2, H-K ~0.0 mag. For a standard extinction law ([Rieke fc Lebofsky 1985] ), 
and the young population colors, the derived visual extinctions are Ay= 17 (W), 15(S), 
24(SE), 18(NE), and 13(A11) mag based on the observed H-K magnitudes given in Table 
1. For an old population (J-H ~0.74, H-K ~0.24 mag), the visual extinctions are 13, 11, 
20, 14, and 9 mag respectively, similar to the 1" resolution estimates obtained by Larkin 
et al. (1995) from Pact and Br7 line ratios. Considerably higher estimates were obtained 
by Sturm et al. (1996) and it is clear that there are undoubtedly areas in the nucleus from 
which virtually no radiation escapes, even at 2.2/im. 
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3.3. Nuclear Structure 



The complex near infrared structure of Arp 220 is due to there being multiple centers of 
star formation activity (as in Arp 299, cf. Wynn- Williams et al. 1991) and strongly varying 
dust obscuration within the merger nuclei. The relatively simple double-nuclei structure of 
the radio continuum (Paan fc Haschick 1995Q and the mm-wave dust continuum QScoville, 



Yun, fc Bryant 1997|) suggest that the major mass concentrations are associated with the W 



and SE sources. The NE and S sources would then be relatively minor mass concentrations, 
albeit with high present-epoch luminosity. 

If the intrinsic, unobscured shapes of the expected stellar emission are either ellipsoidal 
(in the case of the nuclear core clusters), point-like (for unresolved stellar clusters) or 
disk-like (for young stars formed in disk), the unobscured morphology should appear 
spherical or elliptical depending on the viewing angle. For both the E and W nuclei, the 
morphology at 1.1-2.2/zm is clearly not elliptical : the western nucleus appears as a crescent 
and the eastern nucleus is split into two components (NE and SE) with very different colors. 
Given our suggested registration of the radio and infrared, it is evident that the true centers 
of the nuclei are so obscured that neither is seen directly at 2.2/im. 

The crescent or partial-ring morphology of the western component is clearest in the 
deconvolved 2.2/im image (Figure 4), but it is also apparent in the original 2.2/im data 
(Figure 2). This morphology might readily arise under two circumstances : if there is an 
obscuring disk of dust (and gas) embedded within the spheroidal nuclear star cluster or if a 
central starburst ring or disk of young stars is partially obscured by its own dust. Figure 
6 models the surface brightness distribution for a star cluster in which the midplane is cut 
by a totally opaque dust disk; here, the observed crescent-shape arises due to an increase in 
geometrical path length through the cluster periphery. The second hypothesis with a ring 
or torus of massive star formation is similar to models proposed for Seyfert galaxy nuclei 
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(eg. Maolino & Rieke 1995, Genzel est al. 1995). 

Two important characteristics of the observed emission should be emphasized : the 
crescent shape and the extremely abrupt cutoff by dust to the south of the W source. The 
former strongly implies that the dust is embedded in the core cluster/disk rather than being 
a foreground dust lane; the latter implies that the dust is in a extremely thin, but very 
opaque disk. This geometry agrees with that suggested by Scoville, Yun & Bryant (1997) 
based on arcsec resolution imaging of the mm-CO line. The low velocity dispersion of the 
molecular gas implied a very thin disk (Az=16 pc). A similar disk structure in the nuclei 
has been suggested by Baan & Haschick (1995) based on their H 2 CO data. 

The morphology of the eastern nucleus is not so distinctive, yet the apparent bifurcation 
of the nucleus into lowly and highly reddened components (NE and SE respectively) and the 
probable positioning of the radio nucleus in an area of weak 2.2/zm emission both suggest 
a nuclear dust disk - like that in the W nucleus. For both nuclei, the highest extinctions 
occur to the south of the radio centers, suggesting that the dust disks are tilted such that 
the south is the near side. 

The measured infrared sizes (FWHM) of the W, SE and NE source components are 
0.2-0.5" , corresponding to 75-185 pc. The separation of the radio nuclei (A and B) is 0.98" 
or 364 pc in projection. Their actual separation is ~600 pc if the geometry suggested by 
Baan & Haschick (1995) and Scoville, Yun & Bryant (1997) is correct since the nuclei would 
then be well off the major axis of their orbits. Within each of the nuclei, Scoville, Yun & 
Bryant (1997) detected CO emission profiles of full-width ~250 km s~ x . If the rotation 
velocities of gas within each nuclear cluster is ~125 km s _1 at 75 pc radius, their dynamical 
masses are ~3 x 10 8 AL. 
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3.4. Unresolved Clusters 

Magnitudes and positions of eight clusters for which reliable measurements could be 
made are listed in Table 1 in order of decreasing flux at 1.1 /mi. All of these sources 
are too bright to be individual stars. Given our resolution (~ 30 pc), the fact that they 
are unresolved is entirely consistent with their being globular clusters. Their measured 
offsets from the 2.2/zm centroid imply projected galactic radii of 0.9-2.7 kpc. The apparent 
clustering of these sources around Arp 220, with none seen beyond 7.5" radius, suggests 
that most of them must be associated with Arp 220 rather than in the foreground. Given 
the fact that we have measured only two colors and must use one color to constrain the 
reddening, the remaining color does not strongly discriminate the ages of these clusters. 
The cluster sources outside 2" radius have measured colors at 1.1, 1.6, and 2.2/rni which are 
consistent (within the uncertainties) with those of globular clusters in NGC 5128 QFrogel 



1984|) , indicating that they could have age > 10 9 yrs and essentially no near- infrared 



reddening (Ay < 5mag). The absolute K magnitude of the brightest cluster (M K = -13.5 
mag) in Arp 220 is ~1.5 magnitudes brighter than the brightest cluster in NGC 5128 
flKrogcl 1984| ) and is comparable in brightness to the tip of the globular cluster luminosity 
function for ellipticals flAjhar et aL 1994] ). very young much that massive, globular clusters 
event in the Arp 220 (cf. ). 



The NICMOS project has been supported by NASA grant NAG 5-3042. It is a pleasure 
to acknowledge assistance and discussions of L. Armus, L. Bergeron, C. Boone, J. Goldader, 
A. Fruchter, J. Jensen, T. Soifer, J. Surace, and R. White. We are grateful to K. Borne and 
R. Lucas for sharing their I-band image of Arp 220 with us. 
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Fig. 1. — False color image of the 1.1, 1.6, and 2.22/mi data for the entire field of view (19" 
). The long wavelength observations were deconvolved to 0.14" resolution. 

Fig. 2. — Contour and grey-scale images for the 1.1, 1.6, and 2.22/im data for the central 
2.5" are shown together with the ratio 2.22/im/l.l/im using the original-resolution drizzled 
images. The coordinates are offsets in a and 5 from the peak at 2.22/iin (01950= 15 /l 32 m 46 s 90, 
^1950= +23°40'07"94). The contours are spaced logarithmically by factors of 1.16, 1.17, and 
1.16 down from the peak on the W source at 5.1, 15, and 20 mJy (arcsec)~ 2 at 1.1, 1.6 and 
2.22/im. In the ratio map, the peak value is 15.8 and the contours are spaced by factors of 
1.54. 

Fig. 3. — Similar to Figure 2 except that the short wavelength observations were convolved 
to the a resolution of 0.22"( (i.e. the original resolution at 2.22/mi). The contours are spaced 
logarithmically by factors of 1.13, 1.14, and 1.16 down from the peak on the W source at 
3.2, 11, and 19 mJy (arcsec)~ 2 at 1.1, 1.6 and 2.22/im. In the ratio map, the peak value is 
15.2 and the contours are spaced by factors of 1.34. 

Fig. 4. — Similar to Figure 2 except that the long wavelength observations were deconvolved 
to 0.14" resolution using Richardson- Lucy algorithm without damping. The contours are 
spaced logarithmically by factors of 1.21, 1.22, and 1.26 down from the peak on the W 
source. In the ratio map, the contours are spaced by factors of 1.65. 

Fig. 5. — Radio continuum sources at A=6 cm are superposed as contours on the 'extinction- 
corrected' 2.22/im image (see text). The radio source structure is derived from the 
deconvolved source sizes and PA's given by Baan & Haschick (1995). The registration of 
the radio and infrared is discussed in the text. The weak radio component C coincides with 
the infrared source S which is not apparent in the corrected 2.22/im image due to its low 
extinction, ie. relatively blue color. The coordinates are offsets from the peak at 2.22/im 
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(a 1950 = 15 fe 32 m 46 s 90 , 5 1950 = +23°40'07"94). 

Fig. 6. — A galactic-nucleus star cluster with an embedded, opaque dust disk provides a 
possible explanation of the observed morphology in the Arp 220 W source (and possibly the 
NE-SE complex). In the model the stellar density is constant out to the cutoff radius and 
the disk is inclined 20° to the line-of-sight. 
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Table 1. Source Properties 
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a Offscts in " relative to the peak at 2.2^m, ai 950 = 15 h 32 m 46?90, ^1950 = +23°40'07'.'94 based on the registration of 
the infrared and radio as discussed in the text. 

b Magnitudes based on the flux calibration of Rieke et al. (1997) and have typical uncertainty of ±0.1 mag. m = 
corresponds to 1922, 1075, and 667 Jy at 1.1, 1.6, and 2.22^m. 

c Magnitudes interpolated from nii.i and nii. 6 magnitudes. 
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